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SUMMARY 

The  atmospheric  differential  absorption  by  water  vapour 
between  two  C02  laser  transmission  lines  has  been  measured 
up  to  heights  of  1.5  tan  with  a  groundbased  sounder  system. 
Water  vapour  profiles,  some  containing  distinct  layers, 
have  been  estimated  from  long  and  short  pulse  soundings, 
and  compared  with  radiosonde  measurements.  Strongly 
scattering  aerosol  layers,  often  associated  with 
temperature  inversions,  cause  severe  distortion  of  the 
long  pulse  differential  absorption  profiles,  indicating 
the  necessity  of  using  a  short  pulse  technique.  This 
technique  shows  potential  usefulness  for  the  remote 
measurement  of  radio  refractive  index  layers  in  the 
atmosphere. 
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1 .  INTRODUCTION 

Task  DST  77/116  has  as  its  main  aim  the  development  of  a  Carbon  dioxide  (C02) 
infra-red  laser  sounder  to  measure  water  vapour  profiles  in  the  atmosphere. 
This  report  provides  results  of  measurements  up  to  November  1981  and  documents 
the  status  of  the  program  at  that  time.  The  conclusions  presented  will  be  of 
a  tentative  nature,  pending  further  measurements  and  the  publication  of  a 
subsequent  report. 


2 .  BACKGROUND 

Following  the  advent  of  the  laser,  efforts  were  soon  made  to  realise  its 
potential  for  remote  sensing  of  atmospheric  constituents.  Initial  work  was 
mainly  on  measurements  of  the  concentration  and  structure  of  atmospheric 
aerosols  by  the  LIDAR  technique  (Light  Detection  and  Ranging).  Later  this 
technique  was  extended  to  utilise  the  principle  of  Raman  spectroscopy  in  which 
each  molecular  species  reradiates  at  a  frequency  which  is  uniquely  offset  from 
the  exciting  laser  frequency.  However,  because  of  the  low  level  of  the 
reradiated  Raman  energy,  this  technique  is  not  feasible  in  a  daylight 
background. 

Another  technique,  which  utilises  the  frequency  dependent  absorption 
properties  of  water  vapour  on  backscatter  returns  from  atmospheric  aerosols  to 
determine  water  vapour  content  along  the  path,  was  developed  by  Schotland  in 
1964(ref.l).  Although  initially  the  method  appeared  to  be  limited  by 
difficulty  in  achieving  the  necessary  stability  of  the  ruby  laser 
transmission,  these  difficulties  have  now  been  overcome  and  other  high 
stability  lasers  which  are  suitable  for  implementing  the  technique,  such  as 
the  C02  laser,  have  been  developed. 


3.  DIFFERENTIAL  ABSORPTION  AT  LASER  WAVELENGTHS 

A  wavelength  tunable  C02  laser  is  capable  of  tuning  over  many  tens  of  discrete 
lasing  frequencies  which  correspond  to  energy  transitions  in  the  C02  molecule. 

Over  the  C02  band  of  emissions  (9.6  to  10.8  p)  there  is  a  continuous 
background  absorption  due  to  water  vapour.  However,  there  are  also  some 
discrete  absorption  lines  due  to  water  vapour  which  overlap  some  of  the  C02 
lines  and  hence  these  C02  lines  are  more  strongly  absorbed  by  water  vapour. 
The  R(20)  C02  line  (10.2466  p)  is  located  in  a  water  vapour  absorption  line 

and  is  much  more  strongly  absorbed  than  any  other  C02  transmission  line.  This 
line  has  been  used  in  conjunction  with  several  other  lines  having  lower 
absorption  to  measure  the  differential  absorption  along  the  transmission  path. 

The  backscattered  power  detected  by  a  lidar  system  can  be  expressed  as(ref.l). 

Pr  =  Pt(c  fJV  e*p('2  f  (▼  ♦  O) 

Jo 

when 


=  transmitted  power  per  pulse  (W) 
c  =  velocity  of  light(m/s) 


S  =  scintillation  error  due  to  wavelength  variability  of 
scatter. 

In  deriving  equation  (2)  it  has  been  assumed  that  the  differences  in  the 
extinction  and  barkscattering  coefficients  at  the  two  wavelengths  are 
negligible.  This  is  a  reasonable  approximation  for  the  CO2  line  pairs 
employed.  A  careful  selection  of  the  two  wavelengths  used  for  the  absorption 
measurements  will  minimise  the  error,  I,  due  to  absorption  by  other  molecular 
species.  Errors  due  to  scintillation  are  usually  negligible  over  small 
wavelength  ranges. 
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Since  Schot  land's  initial  experiment,  the  technique  has  been  extended  to  the 
near  infra-red  (0.715  to  0.740  p)  with  a  ruby  pumped  dye  laser  and  to  C02 
wavelengths  (10.6  p)  where  suitable  water  vapour  absorption  characteristics 
are  available.  Table  1  is  a  representative  list  of  remote  sampling 
measurements  of  atmospheric  gases. 


4.  SOUNDER  DEVELOPMENT 


4.1  Equipment  description 

A  block  diagram  of  the  sounder  components  is  sketched  in  figure  1  The 
Lumonies  TEA-103-2  laser  together  with  the  30  cm  Newtonian  telescope  is 
supported  on  a  mobile  steel  carriage.  Specif ications  for  the  sounder  are 
given  in  Table  2. 

The  Laboratory  in  which  the  laser  is  housed  has  a  one-metre  square 
ventilation  chute  protruding  through  the  roof.  A  cover  is  attached  to  this 
which  can  be  slid  off  to  one  side  on  rollers  to  open  the  one-metre  square 
aperture.  The  telescope  points  vertically  through  the  centre  of  the 
opening. 

The  C02  laser  is  optically  aligned  approximately  with  the  aid  of  a  helium 
neon  laser,  which  is  inserted  as  indicated  in  figure  1.  More  precise 
alignment  is  subsequently  carried  out  over  a  horizontal  range  with  the  aid 
of  a  target  van.  The  procedure  is  described  in  Appendix  I. 

4.2  Radio  frequency  interference 

Following  laser  and  telescope  alignment,  and  a  demonstration  that  the  laser 
performance  measured  with  a  target  van  was  within  the  specified  performance 
limits,  a  high  sensitivity  mercury-cadmium-tel luride  (CMT)  detector, 
designed  for  operation  at  77°K,  was  substituted  for  the  Pl-61  pyroelectric 
detector.  Severe  radio  frequency  (rf)  pickup  problems  were  then 
encountered,  due  to  the  relatively  large  body  of  the  detector  dewar  and  the 
high-gain,  wide-bandvidth  preamplifier  (5  Hz  to  10  MHz)  which  was  attached 
to  the  dewar. 

It  was  ascertained  that  the  rf  pickup  from  the  35  kV  discharge  of  the  laser 
power  supply  drove  the  preamplifier  into  saturation.  This  was  of  concern 
for  two  reasons.  Firstly,  there  was  subsequent  distortion  to  any  detected 
infra-red  signal  and  secondly,  there  was  was  a  danger  of  damage  to  the 
front  end  of  the  preamplifier,  due  to  excessive  signal  levels.  To  reduce 
the  rf  radiation  from  the  laser,  a  box  shield  was  constructed  around  it,  as 
shown  in  figure  2.  This  effectively  meant  the  laser  was  double  shielded, 
as  its  primary  enclosure  was  specified  as  providing  rf  shielding  by  the 
makers.  Although  this  procedure  reduced  the  rf  pickup  on  the  detector  by 
approximately  20  dB,  there  was  still  an  undesirable  level  of  rf 
breakthrough  present.  To  reduce  this  further,  a  shield  was  constructed 
around  the  detector.  Both  the  laser  and  detector  were  earthed  at  one  point 
only  to  minimise  earth  current  effects.  The  enclosure  for  the  detector  is 
shown  in  figure  3  (the  telescope  is  to  the  right  of  the  detector  shield 
shown).  Copper  braid  was  inserted  between  contacts  on  all  aluminium 
surfaces  to  reduce  radiation  leakage.  The  window  in  the  side  of  the 
detector  housing  is  to  allow  adjustment  of  the  detector  which  is  mounted  on 
a  carriage  with  adjustments  available  along  3  axes. 

4.3  Wavelength  calibration  of  laser 

A  spectrum  analyser  designed  for  use  at  C02  wavelengths  (Optical 
Engineering  C02  Spectrum  Analyser)  was  employed  for  the  calibration.  As 
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the  normal  display  scale  is  not  able  to  withstand  high  power,  it  was 
replaced  by  a  carbon  strip  on  which  the  beam  impinged  adjacent  to  the 
wavelength  scale  markings.  The  scale  calibration  was  initially  checked  by 
using  a  helium  neon  laser  as  described  in  the  handbook.  The  scale  read 
high  by  0.007  microns  and  this  was  allowed  for  in  the  calibration.  A 
special  platform  was  constructed  to  mount  the  spectrum  analyser  and  a 
90  degree  light  tube  was  built  to  direct  the  laser  beam  into  the  aperture 
of  the  spectrum  analyser.  A  carbon  strip  with  a  small  aperture  was  placed 
in  front  of  the  spectrum  analyser  to  protect  areas  adjacent  to  the  aperture 
slit  from  the  CO2  beam. 

The  positions  of  several  C02  wavelength  lines  on  the  C02  laser  grating 
micrometer  scale  were  noted.  It  was  later  found  that  maintenance  work  on 
the  output  optical  system  of  the  laser,  or  severe  vibration  of  the 
micrometer,  alters  the  micrometer  calibration,  so  that  it  is  advisable  to 
recalibrate  in  these  cases. 

4.4  Digital  recording  of  range  intensity  display 

Brief  details  of  all  components  of  the  laser  sounding  and  digital  recording 
system  are  given  in  Appendix  II.  One  of  the  initial  recording  problems 
encountered  was  to  sample  the  lidar  return,  which  has  a  dynamic  range  of 
approximately  40  dB,  at  rates  greater  than  5  MHz.  Because  of  the 
unavailability  of  high  speed  A/D  converters  with  greater  than  6  bit 
sampling  range  at  the  commencement  of  the  programme,  it  was  decided  to 
perform  a  logarithmic  amplification  on  the  signal.  As  the  detector  output 
voltage  is  directly  proportional  to  the  incident  optical  power,  the 
dynamic  range  of  the  detected  signal  (in  dB)  is  twice  that  of  the  optical 
power  (in  dB).  Hence  an  A/D  system  capable  of  measuring  a  40  dB  variation 
in  optical  power  must  be  able  to  operate  over  a  dynamic  range  of  80  dB. 
With  a  6  bit  A/D,  an  electrical  power  resolution  of  1.25  dB  is  obtainable, 
corresponding  to  an  optical  power  resolution  of  0.625  dB.  The  logarithmic 
amplifier  response  deviates  from  a  logarithmic  function  at  the  bottom  and 
top  of  the  ranges  of  operation.  A  calibration  is  shown  in  figure  4, 
indicating  a  true  logarithmic  response  over  a  52  dB  range. 

The  A/D  device  was  also  employed  to  monitor  the  transmitted  laser  pulse, 
and  is  described  in  Appendix  II. 

4.5  Initial  horizontal  transmission  tests  with  CUT  detector 

Following  completion  of  shielding  on  both  laser  and  detector,  the  system 
was  tested  on  the  target  van  after  inserting  40  dB  optical  attenuation  in 
front  of  the  detector.  Strong  echos  were  received  from  the  van  showing  the 
relatively  great  improvement  in  performance  obtainable  from  the  CMT 
detector  compared  to  the  pyroelectric  detector. 

The  echo  strengths  received  from  the  side  of  the  van  were  measured  for 
three  adjacent  C02  lines,  R(22)  (10.233  p),  R(20)  (10.247  p)  and 
R(18)  (10.260  p).  The  measured  water  vapour  absorption  coefficients  for 
R( 18 )  and  R(22)  are  both  low,  and  significantly  less  than  for  R(20)(ref .5) . 
The  water  vapour  content  of  the  atmosphere  near  the  earth's  surface  was 
calculated  from  wet  and  dry  bulb  temperature  observations  taken  near  the 
transmission  position.  The  predicted  and  measured  differential  absorptions 
relative  to  R(20)  were  then  compared.  (The  predicted  differential 
absorptions  for  R(18)  to  R(20)  and  R(22)  to  R(20)  were  very  nearly  equal). 
In  general  the  measured  differential  absorption  for  R(18)  to  R(20)  was 
slightly  lower  than  predicted.  This  may  have  been  due  to  slight  changes  in 
laser  beam  transmission  direction  as  the  wavelength  was  changed  by  tilting 
the  diffraction  grating.  Measurements  of  the  beam  imprint  on  heat 
sensitive  paper  at  10  m  confirmed  this  as  a  likely  explanation.  Another 
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possible,  but  not  likely,  explanation  is  an  additional  differential 
absorption  effect  due  to  another  gas  such  as  ammonia . 

As  the  differential  absorption  was  only  of  the  order  of  2  dB  for  the  above 
measurements  (at  450  m  range),  the  measurements  were  repeated  at  longer 
ranges  using  the  hills  face  as  a  target  (shown  in  figure  5).  Oscilloscope 
photographs  of  echos  are  shown  in  figure  6,  the  range  scale  is  750  m  per 
division.  However,  in  this  case  at  a  range  of  3.5  km  the  absorption  was 
such  that  an  R(20)  echo  could  not  be  detected  above  noise  level  for  single 
pulse  operation.  This  is  to  be  expected,  as  the  differential  absorption, 
in  decibels,  is  directly  proportional  to  range  for  a  constant  water  vapour 
content  along  the  range,  and  therefore  the  absorption  should  be  nearly 
8  times  the  value  for  450  m  range,  ie  a  value  of  15  to  16  dB  would  be 
expected . 

It  was  concluded  that  the  differential  absorption  effect  on  the  sounder 
performance  was  measurable  over  useful  ranges. 


5.  VERTICAL  MEASUREMENTS  OF  DIFFERENTIAL  ABSORPTION 

5.1  Introduction 

During  the  horizontal  transmission  tests  it  was  noted  that  strong 
interfering  echos  were  received  from  the  structure  supporting  the  roof 
mirror.  This  was  due  to  reflections,  by  this  structure,  of  the  forward 
scattered  laser  beam.  These  reflections  were  reduced  by  covering  as  much 
as  possible  of  the  structure  with  a  non-reflecting  material.  However, 
discernible  echos  still  remained.  Finally  to  reduce  echos  to  a  minimum  a 
75  mm  diameter  light  pipe  was  used  to  completely  enclose  the  laser  beam, 
until  it  cleared  the  roof  hatch  as  shown  in  figure  7.  However,  when  the 
roof  mirror  was  employed  for  horizontal  transmissions ,  echos  still 
remained,  probably  due  to  reflections  from  the  roof  hatch,  after  the  beam 
had  been  reflected  by  the  laser  roof  mirror.  It  was  concluded  that  to 
completely  eliminate  these  echos  for  future  horizontal  sounding 
measurements  the  laser  beam  would  need  to  be  completely  enclosed  until  it 
cleared  the  building.  This  step  has  not  been  taken  to  date  because  it  is 
intended  that  the  system  be  mainly  operated  in  the  vertical  mode. 

5.2  Preliminary  measurements 

These  consisted  of  digitally  recording  sequences  of  five-ten  pulses  of  each 
of  the  three  C02  lines  R(18),  R(20)  and  R(22).  The  laser  was  operated  with 
a  gas  mixture  of  He  and  C02  only,  which  gives  the  shortest  pulse  length, 
but  at  a  lower  power  output  than  when  nitrogen  (N2)  is  included  in  the  mix. 
Initial  pulse  echos  from  the  target  van,  with  different  N2  mixes,  indicated 
a  considerable  variation  in  the  pulse  length.  Further  direct  measurea^nts 
on  the  emitted  pulse  showed  a  pulse  length  variation  of  over  10  to  1 ,  with 
the  shortest  pulses  corresponding  to  the  zero  N2  mix.  The  pulse  length  is 
defined  as  the  time  taken  to  decay  to  a  value  of  e-1  of  its  peak  value. 
The  peak  value  of  the  pulse  is  reached  within  50  ns  for  all  gas  mixes,  but 
the  pulse  decay  varies  between  100  ns  for  a  zero  N2  gas  mix,  up  to  1  ps  for 
a  standard  N2  mix. 

5.3  Effects  of  pulse  energy  and  duration 

Preliminary  measurements  indicated  a  range  capability  with  the  short  pulse 
of  800  m  to  1  km  before  the  differential  absorption  rate  began  to 
apparently  decrease  with  range.  This  condition  is  reached  when  the 
differential  absorption  equals  the  signal-to-noise  of  the  signal  suffering 
the  least  absorption,  ie  typically  the  R(18)  transmission.  When  the  long 
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pulse  high  energy  mode  was  employed,  this  range  limitation  was  extended  out 
to  beyond  1.5  km.  However  in  this  case,  measurements  at  ranges  less  than 
500  m  indicated  that  range  gated  samples  were  often  distorted  by  the 
contribution  of  backscattered  energy  in  the  tail  of  the  pulse  at  the 
shorter  ranges.  Pulse  tail  effects  are  considered  in  Appendix  III.  The 
severity  of  this  interference  depends  on  thp  backscattering  aerosol  density 
profile.  It  is  minimised  when  there  is  a  low  aerosol  density  decay  with 
height,  ie  a  large  scale  height. 

5.4  Pulse  averaging  considerations 


Initial  measurements  indicated  there  could  be  a  considerable  difference 
between  successive  sets  of  averaged  differential  absorption  data,  where 
each  set  was  obtained  from  typically  five  pairs  of  pulses.  Normal  random 
variations  were  expected  due  to  the  random  variations  in  aerosol 
distribution,  but  these  differences  indicated  statistically  non-stationary 
processes  were  occurring.  It  was  therefore  decided  to  add  the  capability 
of  switching  back  and  forth  automatically  between  the  selected  pair  of 
differentially  absorbed  wavelengths,  so  that  non-stationary  effects  would 
be  minimised.  In  the  original  concept  of  the  microprocessor  program,  the 
intention  was  that  the  pulse  triggering  and  line  selection  would  be  under 
program  control.  However,  due  to  lack  of  staff  resources,  this  facility 
has  not  yet  been  implemented.  At  present,  periodic  triggering  sequences 
are  generated  by  the  laser  control  unit,  running  at  a  pulse  repetition 
period  of  13  s.  The  drive  for  the  optical  grating  is  also  initiated  by 
this  trigger. 


The  width  of  a  sampled  range  cell  was  25  m,  which  was  therefore  the  maximum 
resolution  of  the  system.  However,  in  processing  the  backscattered 
absorption  data,  it  was  found  to  be  advantageous  to  compute  a  running 
average  of  the  results  over  at  least  ten  25  m  range  cells  to  smooth  out 
quantizing  effects.  Four  different  processing  methods  were  used.  In  the 
first  estimation  method,  the  25  m  range  cell  water  vapour  estimates  in 
grams/kilogram  were  smoothed  over  2n  +  1  range  cells  obtained  from  the 
averaged  power  P  for  N  pulse  pairs  in  adjacent  range  cells,  as  follows: 


Wj 


(1) 


i=j+n 

t 


w. 

i 


i=j-n 


622  j- 

2AKAR 


l0ge(j 


,  Pis(ri) 

Pi»(ri+1)N 

\Pio(r.) 

Pjo(ri+1)/ 

(3) 


where  AK  is  a  constant  determined  by  the 
coefficients  for  each  line  (km-1  atm-1) 
(25  m). 


difference  in  the  absorption 
and  AR  is  the  range  interval 


The  second  processing  method  consisted  of  summing  the  power  ratios  over 
2n+l  cells  before  taking  the  logarithm,  ie: 
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The  third  method  consisted  of  suming  power  estimates  before  ratios  were 
taken,  ie: 
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where 


i-j+n 

£  =  X  ' 

i=j-n 

The  final  method  consisted  of  processing  the  pulse  to  pulse  ratios, 

«j(4)  =2fiW£(Ri  -  Rr+1> 


with  1  having  same  meaning  as  above  and  the  being  summed  over  N  pulse 
pairs,  ie: 

N 

R.  =  \  log  (P18  (r  )/P20  (r.)) 

I/O  1  1 


An  expression  for  the  mean  square  error  of  the  estimates  is  given  in 
Appendix  IV. 

The  averaging  methods  (2),  (3),  and  (4)  are  represented  by  a  cross,  a 
diamond  and  a  triangular  symbol  respectively  in  the  plots  of  smoothed 
differential  attenuation  versus  range  shown  in  subsequent  figures.  Method 
(1)  is  not  considered  in  these  plots,  because  this  provides  only  a 
smoothing  of  the  raw  water  vapour  estimates,  compared  with  the  smoothing  of 
absorption  differentials  provided  by  the  other  three  methods. 

For  the  purposes  of  computing  water  vapour  content,  method  (4)  was 
considered  the  most  reliable  as  it  should  minimise  non-stationary 
transmission  and  scattering  effects  of  the  atmosphere.  This  method  is  the 
one  which  will  be  used  for  comparison  with  radiosonde  estimates. 


6.  ANALYSIS  OF  SOUNDER  RESULTS 

Up  to  the  present  time,  sounder  profiles  have  been  obtained  on  three  separate 
occasions.  These  will  be  covered  in  order  of  significance,  rather  than 
chronologically . 

6.1  Results  for  3/12/81 

6.1.1  Sounder  data 

The  early  morning  radiosonde  flight  from  Adelaide  Airport  (0830  Central 
Standard  Time  (CST))  had  a  12°C  temperature  inversion  at  1  km  with  an 
associated  drop  out  of  moisture  level.  The  water  vapour  content  of  the 
atmosphere  was  of  the  order  of  6  g/kg  up  to  a  height  of  1  km,  then  dried 
out  to  almost  zero.  Profiles  are  given  in  figure  8,  together  with 
associated  profiles  for: 

(a)  refractive  index  (defined  in  Section  7); 

(b)  dry  refractive  index  (mixing  ratio  assumed  zero); 
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(c)  potential  refractive  index  (temperature,  pressure  and  mixing 
ratio  adjusted  for  a  pressure  of  1000  mbars); 

(d)  potential  temperature  (temperature  adjusted  tor  pressure  ot 
1000/mbars) . 

These  conditions  provided  for  the  desired  characteristic  to  test  the 
sounders  ability  to  detect  a  sharp  discontinuity  in  water  vapour 
content.  The  sounder  was  operated  from  0947  to  1017  (CST)  in  the  long 
pulse  high  power  mode  (Run  1)  and  from  1022  to  1054  in  the  short  pulse, 
low  power  mode  (Run  2).  A  sonde  was  released  at  the  sounder  site  at 
0950.  The  reduced  profile  is  shown  in  figure  9.  This  profile  is  very 
similar  to  the  earlier  one  observed  at  Adelaide  Airport. 

The  long  pulse  backscattered  power  as  function  of  range  is  shown  in 
figure  10  for  the  R(18)  and  R(20j  line  pairs.  The  smoothed  differential 
of  the  return  (R(18) -R(20) )  is  shown  in  figure  11  for  the  different 
processing  methods.  The  running  means  were  computed  over  11  adjacent 
points  to  minimise  quantization  effects.  This  will  be  further 
considered  in  Section  6.1.2. 

The  short  pulse  backscattered  power  and  differential  absorption  results 
(figures  12  and  13)  show  a  generally  linearly  increasing  absorption 
difference  up  to  950  m.  This  is  consistent  with  the  fixed  value  of  the 
water  vapour  content  up  to  this  height.  The  absorption  difference  then 
remains  constant  up  to  1050  m,  indicating  a  negligible  water  vapour 
content.  At  greater  heights  the  absorption  difference  is  degraded  by  a 
falling  signal -to-noise  ratio. 

The  long  pulse  differential  absorption  results  (figure  11)  illustrates 
the  degrading  effect  introduced  at  short  ranges  (<500  m)  by  power  from 
the  tail  of  the  pulse  scattered  back  et  ranges  less  than  the  gated 
range.  However,  at  longer  ranges  the  apparent  increase  in  differential 
absorption  is  initially  greater  than  predicted  by  the  model  described  in 
Appendix  III,  and  then  flattens  out  abruptly  at  a  height  of  1100  m. 

Although  the  flattening  of  the  short  pulse  curve  in  figure  13  at  around 
1000  m  is  most  likely  to  be  due  to  the  sharp  reduction  in  water  vapour 
content  at  this  level,  the  same  conclusion  would  not  hold  for  the  long 
pulse  result  in  figure  11  (where  the  flattening  effect  commences  at 
1150  m).  As  shown  in  Appendix  III,  the  long  pulse  (pulse  length  300  m) 
differential  absorption  would  be  expected  to  continue  to  rise  for  some 
distance  above  the  top  of  the  water  vapour  layer,  due  to  the 
contribution  of  the  reflected  power  below  the  top  of  the  layer.  The 
explanation  for  the  failure  of  the  measured  curve  to  follow  this 
prediction  can  be  found  in  figure  10,  which  shows  the  presence  of 
enhanced  backscattering  in  the  region  800  to  900  m,  just  below  the 
inversion.  The  effect  of  this  is  to  increase  the  apparent  gradient  of 
the  differential  absorption  in  the  layer,  and  immediately  above  the 
layer  a  flattening  occurs.  This  is  shown  in  Appendix  IV  by  a  computer 
simulation,  employing  the  model  of  Appendix  III  with  the  addition  of  an 
aerosol  layer. 

6.1.2  Comparison  with  radiosonde  profile 

Figure  14  is  a  replot  of  the  mixing  ratio  profile  previously  shown  in 
figure  9,  compared  with  the  data  from  the  short  and  long  pulse  sounder 
runs  on  3/12/81.  The  short  pulse  sounder  points  show  a  layer 
terminating  at  a  slightly  lower  level  than  the  radiosonde  profile  and 
are  significantly  different  from  that  curve,  allowing  for  errors  in  each 
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estimate  (rms  error  bounds  are  plotted).  Profile  points  for  the  long 
pulse  sounding  are  not  shown  above  1100  m  because  they  were  considered 
to  be  too  seriously  affected  by  the  aerosol  layer  previously  mentioned. 
They  have  also  been  omitted  at  levels  below  700  m  because  of  the 
severity  of  pulse  tail  effects  in  this  region.  To  investigate  the 
source  of  the  discrepancy  between  the  short  pulse  and  radiosonde 
profiles  a  number  of  unsmoothed  short  pulse  differential  absorption 
profiles  were  inspected.  The  profiles  all  exhibited  the  same 
oscillatory  character.  This  is  evident  in  figure  15,  where  averaged 
differential  absorption  values  for  each  range  cell  are  plotted,  based  on 
processing  method  (4),  together  with  a  smoothed  curve  over  11  points. 

The  initial  ascending  section  of  the  unsmoothed  differential  absorption 
curve  in  figure  15  exhibits  the  undesired  effects  of  quantization  by  the 
A/D  converter,  manifested  by  a  staircase  appearance  for  the  trace. 
Unfortunately,  this  limits  the  range  resolution  obtainable  for  the  water 
vapour  estimates,  since  they  are  given  by  the  gradients  of  the 
differential  absorption. 

6.2  Results  for  3/11/81 

6.2.1  Sounder  data 

The  morning  radiosonde  flight  at  Adelaide  Airport  had  an  adiabatic  lapse 
rate  up  to  900  m  (see  figure  16),  with  a  stable  layer  above  this  level. 
The  water  vapour  content  was  4  to  5  g  up  to  800  m,  with  a  dip  at  1450  m, 
then  started  increasing  again  above  this  level. 

The  lidar  was  operated  in  long  pulse  and  short  pulse  modes  between  1016 
and  1102  CST.  The  averaged  backscattered  power  levels  returned  for  the 
R(18),  R(20)  line  pair  in  the  short  pulse  and  long  pulse  modes  are  shown 
in  figures  17  and  18.  The  latter  diagram  shows  the  presence  of  cloud  at 
3400  m  altitude  with  a  12.5  dB  differential  between  the  two  lines  at 
this  height.  The  differential  absorption  results  are  plotted  in 
figures  19  and  20.  A  switching  transient  in  the  input  to  the  A/D 
converter  is  evident  in  both  plots,  in  the  range  interval 
3750  to  4250  m.  This  transient  is  due  to  the  A/D  converter  being 
switched  to  the  monitor  detector,  to  record  the  level  of  the  laser 
pulse . 

The  rapidly  climbing  section  of  the  long  pulse  curve  (figure  20)  is 
again  due  to  a  reduction  in  the  decay  rate  of  the  aerosol  density  at 
about  700  m,  which  is  evident  from  the  backscattered  power  versus  range 
in  figure  18.  However,  this  layer  has  a  much  greater  thickness 
(extending  to  approximately  1400  m)  than  was  exhibited  by  the  aerosol 
layer  present  in  figure  10. 

6.2.2  Comparison  with  radiosonde  profile 

The  profile  comparison  is  made  in  figure  21.  The  long  pulse  profile  is 
in  reasonable  agreement  above  800  m,  with  lower  altitude  estimates 
seriously  affected  by  the  long  pulse  tail  effect.  The  spuriously  high 
value  of  the  mixing  ratio  at  700  m  is  probably  due  to  the  appearance  of 
an  aerosol  layer  at  about  this  height.  The  error  bars  in  figure  21  are 
generally  greater  than  for  figure  14  because  the  backscattered  power 
showed  a  wider  fluctuation  on  this  occasion.  The  smoothing  procedure 
employed  on  all  profiles  has  removed  much  of  the  oscillatory  cosq>onent 
(thought  to  be  due  to  quantization  noise),  which  is  also  particularly 
noticeable  in  other  short  pulse  profiles. 


RL-0230-TM 


10 


6.3  Results  for  24/11/81 

6.3.1  Sounder  data 

Three  sounder  runs  were  performed  on  this  day,  runs  one  and  two  using 
the  long  pulse  mode  commencing  at  1002  and  1033  CST  respectively,  and  a 
short  pulse  run  starting  at  1112  CST.  The  DR CS  radiosonde  profile  was 
obtained  at  1112  CST.  Since  the  two  long  pulse  runs,  consisting  of  54 
and  61  pulse  pairs  respectively,  gave  very  similar  results  it  was 
decided  to  combine  the  two  runs  together  to  form  a  single  data  set.  The 
backscattered  power  and  differential  absorption  plots  for  the  long  pulse 
run  are  shown  in  figures  22  and  23  respectively.  The  absorption 
difference  variation  with  range  shows  three  characteristic  regions.  The 
initial  low  gradient  section  is  again  due  to  the  pulse  tail  effect.  The 
rapidly  climbing  section  is  due  to  a  reduction  in  the  decay  rate  of  the 
aerosol  density  at  about  800  m.  This  layer  extends  up  to 
1300  to  1400  m.  The  evidence  for  this  is  found  in  the  average 
backscattered  power  versus  range  results  of  figure  22. 

6.3.2  Comparison  with  radiosonde  profile 

The  resulting  profile  in  the  region  600  m  to  1300  m  is  shown  in 
figure  24,  together  with  the  corresponding  radiosonde  profile.  The 
higher  values  of  mixing  ratio,  relative  to  radiosonde  values,  obtained 
between  700  to  900  m,  are  again  due  to  the  influence  of  the  aerosol 
layer  mentioned  above. 

The  short  pulse  results  are  shown  in  figures  25  and  26.  The 
differential  absorption  measurement  (figure  26)  is  limited  by  signal-to- 
noise  restrictions  above  750  in.  The  oscillatory  component,  arising  from 
quantization  noise,  has  been  removed  from  the  lower  altitude  section  of 
the  short  pulse  profile  of  figure  26  by  smoothing,  as  for  the  two 
previous  cases  considered  in  Sections  6.1.1  and  6.1.2. 

6.4  Effect  of  signal-to-noise  degradation  on  layer  determination 

It  is  apparent  in  the  plots  of  differential  absorption  versus  range,  that 
all  curves  monotonically  increase  at  lower  heights,  when  there  is  good  S/N 
ratio,  and  eventually  begin  descending  as  the  S/N  ratio  falls  below  the 
absorption  differential.  Care,  therefore,  must  be  exercised  to  avoid 
interpreting  a  reduction  in  S/N  ratio  as  a  decrease  in  water  vapour 
content.  Water  vapour  estimates  were  only  calculated  in  the  plots  of 
figures  14,  21  and  24  when  the  signal-to-noise  of  the  weaker  return  was 
above  a  value  of  approximately  10  dB.  Using  power  samples  with  S/N  ratios 
below  this  value  will  lead  to  serious  errors. 

An  illustration  of  this  point  is  provided  by  the  differential  absorption 
curve  of  figure  26  which  has  a  plateau  above  750  m  altitude.  As  previously 
discussed,  this  is  usually  indicative  of  a  very  low  water  vapour  content. 
However,  in  this  case  inspection  of  the  backscattered  power  return 
(figure  25)  shows  that  at  the  plateau  base  height  (6750  m)  the  S/N  ratio  of 
the  R(20)  return  was  less  than  10  dB,  with  evidence  of  corruption  by  noise. 
Hence  the  differential  absorption  curve  was  not  used  for  water  vapour 
content  estimation  in  this  region. 
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7.  FUTURE  DEVELOPMENTS 

7.1  Increase  in  sampling  resolution 

As  discussed  in  Section  6.1.1  the  resolution  obtainable  for  water  vapour 
estimates  was  limited  by  the  six  bit  A/D  converter.  This  converter  is 
presently  being  replaced  by  an  eight  bit  device  and  it  is  anticipated  that 
this  will  significantly  improve  the  resolution  and  accuracy  of  the 
estimates . 

7.2  Improvements  to  signal-to-noise  performance 

7.2.1  Heterodyne  detection 

Heterodyne  detection  improves  the  signal  detection  threshold  by  several 
orders  of  magnitude  over  direct  detection.  This  can  be  implemented  by 
using  a  separate  local  oscillator  which  beats  with  the  backscattered 
energy  and  the  transmitter  to  produce  two  rf  signals  whose  difference 
can  be  controlled  by  a  loop.  The  difference  signal  is  then  envelope 
detected  to  give  the  amplitude(ref . 6) .  A  further  possible  advantage  to 
be  gained  is  the  rejection  of  the  tail  of  the  laser  pulse  by  the  loop 
filter  since  the  decaying  part  of  the  pulse  exhibits  a  changing 
frequency  shift. 

7.2.2  Background  limited  detection 

At  present  the  limiting  noise  in  the  detection  system  arises  from  the 
front  end  of  the  preamplifier  (noise  figure  4  dB) .  It  may  be  possible 
to  reduce  the  amplifier  noise  level  so  that  the  background  sky  noise 
becomes  the  limiting  noise.  Further  noise  reduction  would  then  be 
possible  by  using  an  optical  filter. 


8.  CONCLUSIONS 

The  infra-red  tunable  sounder  is  capable  of  measuring  water  vapour  profiles  to 
at  least  1.5  km,  but  the  height  capability  is  dependent  on  the  atmospheric 
aerosol  distribution.  In  cases  of  strong  temperature  inversions,  the  aerosol 
density  often  builds  up  below  the  inversion  and  then  falls  off  rapidly  above 
the  inversion.  This  can  lead  to  serious  distortion  of  the  differential 
absorption  estimates  for  long  pulse  transmissions.  However,  for  short  pulse 
emissions  the  range  capability  of  the  lidar  may  be  improved  (as  in  the  case  of 
the  measurements  on  3/12/81  due  to  the  increase  in  the  intensity  of  the  return 
from  below  the  inversion  level. 

The  results  in  this  paper  show  the  necessity  of  using  a  short  pulse  emission 
in  the  presence  of  atmospheric  inversions  to  enable  the  features  of 
atmospheric  water  vapour  layers  to  be  resolved.  The  data  presented  show  that 
the  sounder  is  capable  of  profiling  in  the  height  range  where  radio  ducts  are 
most  f requent(ref . 7) ,  and  most  seriously  affect  surface  based  radio  frequency 
sensors  and  communications  systems.  Further  development  and  refinement  of  the 
sounder  S/N  performance  is  required  to  utilize  to  the  full  its  undoubted 
potential,  which  could  rival  the  radiosonde  in  accuracy. 
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TABLE  1 .  RESULTS  OF  DIFFERENTIAL  ABSORPTION  MEASUREMENTS 


Reference 

Molecule 

Wavelength 

(pm) 

Path 

Direction 

Resolution 

(m) 

*  Bowell  et  al  (1979) 
(ref. 2)  NASA 

H20 

0.7243 

Vertical 

100 

Asai , Itabe, 

Igarashi  (1979) 

(ref. 3)  RRL  Japan 

o3 

10.6 

Horizontal 

500 

Murray  et  al  (1976) 
(ref. 4)  SRI 

H20 

10.6 

Horizontal 

30 

Schotland  (1964) 

(ref . 1 ) 

h2o 

0.6943 

Vertical 

300 

'•Simul  taneous  transmission  on  two  laser  wavelengths. 
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TABLE  2.  C02  SOUNDER  SYSTEM  PARAMETERS 


C02  TEA  Laser 

Transmitted  Energy 

(a)  Maximum 

J 

(b)  Without  N2  in  gas  mix 

0.  1  J 

Lines 

R(  18) ,  R(20) ,  R(22)  on  10  pm  band 

Pulse  Width 

(a) 

1  ps  (to  e-1  level) 

(b) 

100  ns  (to  e-1  evel) 

Beam  Divergence 

1  x  0.6  mrad 

Typical  PRF 

0.07  Hz 

Maximum  PRF 

0.5  Hz 

Receiver 

Telescope  diameter 

30  cm 

Field  of  viev 

1  mrad  FWHM 

Hg-O.d-Te  photoconduct ive 

detector 

In 

D* 

3  x  10  cm  Hz  /W 

Time  Constant 

50  ns 

Size 


1  x  1 


0 


Figure  1.  Block  diagram 


ot  LIDAR  and  recording  system 


[iMtli] 


Figure  3.  Double  shielding  enclosure  for  detector 
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PRINT  RECORD  READING  (HEXADECIMAL) 


0  10  20  30  AO  50  60  70 


PRINT  RECORD  READING(DEOMAL) 


Figure  4.  Laser  signal  processing  calibration 
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(a) 


(b) 


Figure  6.  Echoes  from  hills  face  at  three  different  C02  transmission  lines 
(a)  R(18),  (b)  R(22) ,  (c)  R(20) 


HEIGHT  IN  METRES 


HE l &H T  IN  METRES 


•  R  (20)  DATA 
+  R  (18)  DATA 

NUMBER  OP  PULSE  PAIRS  =  65 


* 
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Figure  10.  Returned  power  versus  range  for  long  pulse  run  (at  0947  CST)  on 
3-12-81,  for  R(20)  and  R(18) 
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KK 
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Figure  H.  Differential  absorption  versus  range  for  long  pulse  run  3-12-81 
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Figure  12.  Returned  power  versus  range  for  short  pulse  run  (at  1022  CST)  on 
3- 12-81 ,  for  R(20)  and  R(18) 
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IS-  Average  25  ■  range  cell  differential  absorption  data  for  short 
pulse  run  on  3-12-81 
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Figure  16.  Adelaide  Airport  radiosonde  profiles  0830  CST,  3-11-81 
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Figure  18.  Returned  power  versus  range  for  long  pulse  run  (at  1108  CST)  on 
3-11-81  for  R(20)  and  R(18) 
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RANGE  (METRES) 

Figure  20.  Differential  absorption  versus  range  for  long  pulse  run  on  3-11-81 
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Figure  22.  Returned  power  versus  range  for  long  pulse  run  (at  1002  CST)  on 
24-11-81,  for  R(20)  and  R(18) 
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Figure  23.  Differential  absorption  versus  range  for  long  pulse  run  on  24-11-81 
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Figure  25.  Returned  power  versus  range  for  short  pulse  run  (st  1112  CST)  on 
24-11-81,  for  R(20)  snd  R(18) 
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APPENDIX  I 

OPTICAL  ALIGNMENT  PROCEDURES 

1. 1  Approximate  method 

The  Newtonian  telescope  is  designed  to  image  on  a  detector  mounted  on  the 
side  of  the  telescope.  The  assembly  holding  the  45°  telescope  mirror  B  and 
a  locussing  lens,  also  houses  the  45°  mirror  (hereafter  called  laser  mirror 
A)  which  deflects  the  beam  from  its  initial  horizontal  path  to  a  vertical 
t  ra jectory . 

A  helium  neon  (HeNe)  laser  is  substituted  for  the  C02  laser  by  interposing 
a  mirror  C  at  45°  along  the  patch  between  the  laser  and  laser  mirror  A. 
The  HeNe  beam  is  brought  in  from  the  side,  as  indicated  by  the  dashed  line 
in  the  figure. 

The  Newtonian  telescope  is  aligned  relative  to  the  reference  HeNe  laser 
beam  by  the  following  technique: 

(i)  A  reference  transmission  beam  line  is  established  (using  the  HeNe 
laser),  which  after  reflection  from  mirror  A  passes  through  the  centre 
of  the  roof  aperture.  The  reference  beam  is  centred  on  a  chosen 
position  in  the  plane  of  the  roof  opening  by  adjusting  laser  mirror  A. 
This  is  the  only  adjustment  for  which  this  mirror  is  employed. 

(ii)  The  telescope  is  then  aligned  to  this  reference  by  using  a 
ret rorel 1 ect or  to  reflect  the  reference  ray  back  into  the  telescope  and 
rotating  the  ret roref lector  through  360°  to  generate  a  cylinder  of 
returning  rays  (in  temporal  sequence).  The  axis  about  which  the 
cylinder  of  rays  is  focussed  is  adjusted,  by  means  of  mirror  B,  to  pass 
through  the  centre  of  the  detector.  (The  HeNe  focussing  image  is 
displaced  171  mm  beyond  the  position  of  the  C02  focussing  plane  because 
of  the  wavelength  dispersion  characteristic  of  the  focussing  lens.) 

(iii)  The  HeNe  laser  is  then  removed  and  C02  laser  brought  into 
alignment  with  the  reference  line. 

This  is  not  a  precise  procedure  because  it  relies  on  centering  the  C02 
laser  beam  on  the  260  mm  length  light  tube  (45  mm  diameter)  leading  to 

laser  mirror  A,  and  also  locating  the  beam  on  the  reference  position  on 

the  roof.  However,  it  serves  as  a  starting  point  for  a  more  accurate 
a  1 ignment . 

1.2  Optical  alignment  method 

Once  the  laser  and  telescope  are  aligned  as  in  the  approximate  method,  a 
more  accurate  alignment  can  then  be  achieved  by  the  following  procedure: 

(i)  A  large  mirror  is  mounted  on  the  roof  in  the  centre  of  the 
aperture,  at  45°  to  the  vertical,  to  deflect  the  field  of  view  and 
transmission  back  to  the  horizontal  direction.  This  mirror,  shown  in 
figure  AI-1,  is  constructed  by  attaching  a  75  mm  gold  plated  laser 
mirror  to  a  300  x  180  mm  aluminium  back  plate  via  a  laser  mirror  holder 
which  can  be  adjusted  relative  to  the  back  plate.  An  aluminium  coated 

glass  mirror  with  a  90  mm  hole  in  its  centre  is  then  attached  to  the 

same  back  plate  using  small  blobs  of  silastic  between  the  two  surfaces, 
to  allow  for  differential  thermal  expansion.  The  plane  of  the  central 
laser  mirror  (hereafter  called  laser  roof  mirror)  is  then  adjusted  so 
that  its  plane  is  parallel  to  the  roof  mounted  large  rectangular  glass 
mirror  used  for  collecting  light  for  the  telescope.  Although  this  glass 
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mirror  deviates  from  a  plane  surface  by  several  wavelengths  across  its 
surface  at  the  C02  frequency,  this  will  not  affect  the  alignment 
procedure . 

(ii)  Once  the  roof  mirror  system  is  set  up,  the  laser  beam  is  directed 
at  a  target  van  located  450  m  away.  The  laser  beam  is  centered  on  the 
target  van  by  using  a  small  75  mm  loudspeaker  to  sense  the  laser  pulse 
impact.  The  light  energy  is  converted  to  sound  energy  by  thermal 
expansion  of  the  surface  on  which  the  beam  impacts.  The  laser  beam 
section  is  mapped  out  approximately  by  using  the  acoustic  output  pulse 
level  of  the  speaker  as  a  measure  of  beam  intensity.  At  the  target  van 
the  beam  shape  will  be  found  to  be  elongated  in  the  horizontal  direction 
with  two  hot  spots  horizontally  displaced.  This  is  in  accord  with  the 
predicted  diffraction-limited  beam  section  which  has  a  width  1.7  times 
the  height.  The  telescope  is  then  visually  centered  on  the  target  van 
reference,  using  a  perspex  disc  with  a  crosshair  located  at  the  detector 
position  (axially  adjusted  to  allow  for  wavelength  difference). 

(iii)  After  the  initial  alignment  (done  at  visible  wavelengths  with  the 
telescope)  a  1  x  1  iran  pyroelectric  detector  (Molection  Pl-61)  is  placed 
at  the  C02  laser  far  field  focussing  position.  This  is  done  to  avoid 
the  possibility  of  damage  to  the  mercury-cadmium-telluride  (CMT) 
detector  by  exposure  to  high  power  infra-red  reflections  during  the 
alignment,  and  the  necessity  of  cooling  the  detector  to  77°K,  which 
involves  filling  the  dewar  flask  with  liquid  nitrogen.  Calculations 
indicate  that,  using  a  pyroelectric  detector  which  has  a  detectivity  of 
5  x  I09  cn h/Hz/W  and  bandwidth  of  approximately  100  Hz,  an  echo  should  be 
received.  Following  the  initial  laser  alignment,  an  echo  could  not  be 
detected  at  first  with  this  procedure,  so  a  metal  corner  reflector  was 
constructed  with  side  length  750  mm  to  facilitate  alignment.  This  was 
placed  in  the  centre  of  the  target  van  and  strong  echos  were  received. 
Later,  when  the  alignment  procedure  was  repeated,  it  was  demonstrated 
that  discernible  echos  could  be  received  directly  from  the  side  of  the 
target  van  (a  matt  white  painted  metal  surface),  in  agreement  with  the 
initial  prediction. 
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APPENDIX  II 

ELECTRONIC  CONTROL,  SIGNAL  PROCESSING  AND  DATA  LOGGING  SYSTEMS  (REFER  TO 

FIGURE  1 1 - 1 ) 


11. 1  Laser  control  unit 

(a)  Provides  electrical  power  supplies,  trigger  pulse  and  PRF  control 
to  CO2  laser. 

(b)  Mixes  CO2,  Helium  and  Nitrogen  gases,  and  regulates  flow  to  laser. 
Also  supplies  nitrogen  gas  to  laser  HV  trigger  spark  gap. 

11. 2  C02  Laser 
Consists  of: 

(a)  High  voltage  trigger  power  supply  and  spark  gap. 

(b)  Three  section  C02  laser  with  energy  storage  capacitors. 

(c)  Front  and  rear  optics  with  micrometer  adjustable  wavelength  grating 
at  rear. 

1 1. 3  Automatic  laser  tuning  system 

Automatically  drives  laser  tuning  grating  alternately  between  any  two 
selected  wavelengths  after  each  transmitted  pulse.  The  system  consists  of: 

(a)  A  programmable  counter  unit  with  internal  clock  source.  After  each 
transmit  pulse  this  unit  provides  a  pre-selected  number  of  output  pulses 
and  a  direction  reversal  signal  for  stepping  motor  control. 

(b)  A  stepping  motor  drive  unit  which  provides  power  supplies  and  logic 
for  stepping  motor  drive. 

(c)  Stepping  motor  and  gear  assembly  which  is  mechanically  coupled  to 
laser  tuning  micrometer. 

11.4  Remote  tuning  position  indicator 

Consists  ot : 

(a)  A  photo-coupler  and  amplifier  mounted  adjacent  to  the  laser 
wavelength  tuning  micrometer  drive  gear. 

(b)  A  remote  LED  indicator  which  is  illuminated  when  a  reflective  strip 

on  the  tuning  drive  gear  aligns  with  the  photo-coupler.  This 

corresponds  to  the  reference  wavelengths. 

11.5  Monitor  detector 

Samples/detects  laser  output  pulse  and  provides  an  amplitude  reference 
level  for  each  transmitted  pulse. 

11. 6  Signal  detector 


(a)  Receives/detects  signal  backscatter  from  atmosphere  via  telescope. 
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(b)  The  logarithmic  amplifier  is  employed  to  effectively  increase 
dynamic  range  by  compressing  signal  amplitude  to  match  input  range  of 
A  to  D  converter. 

11.7  Signal  gate,  switch  and  sum 

Performs  gating  of  monitor  (transmission  sample)  and  received  backscatter 
signals  and  provides  a  sequential  (time  multiplexed)  output  of  same  via  a 
summing  and  amplifier. 

11.8  Fast  A  to  D  converter  and  random  access  memory 

(a)  A  6  bit  parallel  A  to  D  converter  converts  combined  monitor  and 
received  signals  to  a  digital  format. 

(b)  The  converted  signal  is  stored  in  random  access  memory. 

(c)  The  6  MHz  clock  generator  determines  sampling  rate  and  range 
resolution  of  25  m. 

11.9  Microprocessor  system 

(a)  The  microprocessor  system  is  based  on  an  "Intel"  8080  CPU. 

(b)  The  microprocessor  is  used  to  control  and  interface  the  operation 
of  peripheral  units  and  to  process  and  output  the  collected  data. 

(c)  A  video  display  unit  and  keyboard  allow  some  interaction  between 
the  operator  and  the  data  collection  system. 

(d)  The  programme  is  contained  in  erasable  read  only  memory,  providing 
ease  of  modification  and  change. 

11.10  Data  recording 

Collected  data  comprises  260  samples  per  pulse  (signal  followed  by  monitor) 
and  time  of  day,  and  is  recorded  as  follows: 

(a)  On  punch  paper  tape  (ASCII  code) 

(b)  On  printer  (hexadecimal) 

(c)  On  storage  oscilloscope  and  photographic  film  (analog). 

NOTE:  Reference  8  provides  a  complete  technical  description 
of  the  sounder  and  data  logging  systems  with  relevant 
drawings  and  circuit  diagrams  included. 
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APPENDIX  III 
PULSE  TAIL  EFFECTS 

If  the  pulse  is  considered  to  comprise  a  sharp  leading  edge  followed  by  an 
exponential  decay  of  the  form  exp(*t/TAU),  where  t  is  the  time  delay  after  the 
leading  edge  and  TAU  is  the  pulse  length,  then  when  the  leading  edge  is 
returned  from  a  range  R,  a  point-on. the  trailing  edge  of  the  pulse  at  x, 

corresponding  to  a  time  interval  — =£ - ,  will  arrive  back  at  the  receiver  at 

the  same  time  as  the  leading  edge.  The  total  energy  received  at  a  given 
instant  from  all  backscattered  components  of  the  pulse  is  proportional  to 


i  exp  - 

Jx=0 


2  (R-x)  ; 
TAU  x  Cj 


/  x2  dx 


where  HS  is  the  scale  height  of  the  aerosol  density,  W  is  the  water  vapour 
concentration  (mixing  ratio)  profile  and  Of  is  the  absorption  coefficient  of 
water  vapour  on  the  absorbed  wavelength  emission.  Combining  this  expression 
with  one  for  the  unabsorbed  transmission  yields  the  differential  absorption 
ratio 


2 (R-x)  1 
TAU  x  Cj 

expj 

r  *  1 

L  HSJ 

i  t 

r  exp' 

1  *i 

f  2a  XW  j 
•  0.622  1 

|  /x2dx 

2  (R-x)  1 

TAU  x  C  1 

> 

1  1 

>  exp 
1  “ 

r  -> 

X  | 

'  HS  J 

/x*dx 

This  expression  has  been  evaluated  numerically  for  several  values  of  TAU  and 
HS.  The  water  vapour  mixing  ration,  W,  was  assumed  to  be  6  g/kg  up  to  a 
height  of  1  km  and  zero  at  greater  heights.  The  curves  in  figure  III.l  show 
the  reduction  in  the  returned  power  differential  (in  decibels)  at  short  ranges 
as  the  pulse  length  is  increased.  Pulse  lengths  of  30  m  and  300  m  are 
representative  of  the  C02  laser's  short  and  long  pulses  respectively. 


Figure  III.l  illustrates  that  as  the  pulse  length  is  increased  the  height  at 
which  the  full  differential  absorption  (due  to  the  water  vapour  layer)  is 
reached,  suffers  progressive  delay. 


Figure  III.l  Predicted  differential  absorption  versus  range  for  pulse  lengths 
representative  of  C02  laser  output 
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APPENDIX  IV 

COMPUTATION  OF  RMS  ERRORS  IN  WATER  VAPOUR 

The  expression  for  the  average  mixing  ratio  (in  grams/kilogram)  between  the 
ith  and  (i+l)th  range  cells  may  be  written 


W(i) 


-622_  P_i  (i  +  l)P;  (i) 

2AKAR  l\pt  (i)  Pj  (i+l)j 


where  Ar  =  R^+^  -  R;  P, (i)  is  the  average  power  received  (over  N  pulses)  for 
a  transmission  wavelength  A!  and  (km-1 atm-1 ) . 

The  error  for  W(i)  may  be  written 


dw  dw 

Aw^  -  apTci)  Ap,(i)  +  Ap,(i+1) 


622  (  Ap2  (j)  AP,  (i-H)  Ap 2  (j+1)  Ap,  (j)  , 
2AkAR  (_  Pi  (i)  Pi  (i+1)  ‘  P2(i+1)  Pi  (i)  J 


Therefore  the  rms  error  in  W  is  given  by 


Aw 

rms 


622  j  Ap, 2  (j)  AP;2  (j+l)  AP, 1  (i-i-1)  APz2  (i)~,  ** 
2ARAK Pi  (i)  +  Pj2  (i+1)  +  Pi1  (i+1)  +  P i*  (i) 


For  computational  purposes  the  mean  square  error  over  2n+l  range  cell 
estimates  can  be  written 


(Aw.)  2 

J 


_  6222  '  a>,K2m 

4Ar5Ax5  (2n+ 1)  ^(2Pi  s  (i) ) 


j?  10-lCiJi, 
(SP»o  (i)V 


♦ 


2oi  g*(i*l) 

(ZP.»(i*l)) 


T 


♦ 


(SPio  (i+l) )5  J 


where  a2  is  the  variance,  with  the  mean  removed,  of  the  average  power 
estimates  in  the  25  a  range  cells,  and  I  is  as  defined  for  equations  (1),  (2), 
(3)  and  (4). 
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THE  EFFECT  OF  AN  ELEVATED  AEROSOL  LAYER  ON  THE  MEASUREMENT  OF  A  DIFFERENTIAL 

ABSORPTION  PROFILE 

The  differential  absorption  model  introduced  in  Appendix  III  has  been  employed 
to  predict  the  effects  of  an  elevated  aerosol  layer  on  differential  absorption 
measurements.  The  curve  in  figure  V.l  was  obtained  using  the  following 
parameters : 

Normalized  aerosol  backscatter  function,  0  =  exp(-h/HS) 

Layer  backscatter  function,  0  =  30  exp(  -jq^) 

Layer  position  =  700  to  1000  m 

Background  scale  height,  HS  =  200  m 

Water  vapour  mixing  ratio  profile  W  =  6  g/kg  0<  h  <  1000  m 

=0  h  >  1000  m. 

The  aerosol  backscatter  function  is  also  shown  in  figure  V.l.  In  addition  to 
the  above  parameters  the  inverse  range  squared  dependence  of  the  integrands  of 
the  integrals  in  Appendix  III  were  assumed  to  commence  at  a  range  of  50  m.  In 
the  range  interval  0  to  50  m  the  inverse  range  squared  function  was  replaced 
by  a  linearly  increasing  function.  This  allows  for  masking  effects  of  the 
telescope  in  this  region  and  also  closely  predicts  the  actual  returned  power 
distribution  as  a  function  of  range. 

The  differential  absorption  profile  which  was  obtained  has  three  main 
features.  The  initial  slowly  changing  section  is  mainly  determined  by  the 
combination  of  the  pulse  length  and  scale  height  parameters  and  hence  the 
gradient  of  this  section  can  accordingly  be  varied  by  changing  these 
parameters.  Reducing  HS  reduces  the  differential  absorption  at  all  levels. 
Secondly,  the  increasing  section  commencing  at  700  m  is  due  to  the  influence 
of  the  aerosol  layer  which  gives  rise  to  an  enhanced  contribution  from  the 
leading  edge  of  the  pulse.  Finally,  the  flattening  of  the  curve  occurs  when 
the  leading  edge  of  the  pulse  emerges  from  the  top  of  the  layer  and  the  return 
from  the  tail  of  the  pulse  is  enhanced. 

The  dashed  curve  in  figure  V.l  is  for  a  short  pulse  sounding  and  indicates  the 
error  contained  in  the  long  pulse  result.  There  is  only  a  slight  kink  in  this 
curve  at  700  m,  at  the  base  of  the  elevated  aerosol  layer. 
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